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Abstract 

We discuss chiral symmetry breaking critical points from the perspective of PCAC, correlation length 
scaling and the chiral equation of state. A scaling theory for the ratio Rt, of the pion to sigma masses 
is presented. The Goldstone character of the pion and properties of the longitudinal and transverse chiral 
susceptibilities determine the ratio i?^ which can be used to locate critical points and measure critical indices 
such as (5. We show how PCAC and correlation length scaling determine the pion mass' dependence on the 
chiral condensate and lead to a practical method to measure the anomalous dimension rj. These tools are 
proving useful in studies of the chiral transition in lattice QED and the quark-gluon plasma transition in 
lattice QCD. 
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1. INTRODUCTION 



Chiral symmetry phase transitions are particularly interesting because of their physical relevance and 
their theoretical constraints. The order parameter for the high temperature quark-gluon transition of QCD 
is the chiral condensate and considerable effort both in phenomonological studies and computer simulations 
are underway [1]. Non-compact lattice QED experiences a chiral transition at strong coupling and one 
needs new practical theoretical tools to understand the character of this transition [2] . Four Fermi models 
with N species of fermions in less than four dimensions have chiral transitions, and 1 /N expansions indicate 
that these transitions correspond to ultra-violet fixed points of the renormalization group [3], [4]. These 
models illustrate the intimate relation between renormalizability, compositeness and hyperscaling in four 
Fermi and Yukawa models. Composite models of the Higgs mechanism which underlies the Standard Model 
typically rely on a chiral symmetry breaking phase transition to make composite mesons which can produce 
heavy W and Z particles. Although the original versions of this idea (Technicolor) copied the asymptotically 
free dynamics of QCD, other versions imagine theories which arc strongly coupled at short distances [5]. 
Interesting models of a heavy top quark based on gauged Nambu Jona Lasinio models have been considered 
[6]. In this paper we shall be particularly interested in models where chiral symmetry breaking is a short 
distance phenomenon and we shall concentrate on non-asymptotically free dynamics. 

Chiral symmetry breaking phase transitions have important simplifying features which will allow us to 
study them in several complementary ways. Our theoretical framework will be that of statistical mechanics 
and concepts borrowed from ferromagnetic phase transitions such as the equation of state, scaling variables, 
correlation length scaling and hyperscaling relations between critical indices, etc. will be used heavily here. 
This 'standard approach' gains far greater predictive power when current algebra relations following from 
the underlying continuous chiral symmetry of the system's Lagrangian arc developed. The most familiar 
deduction of this approach is the Goldstone character and low energy theorems concerning the pion. But 
we shall present additional equally firm results here involving the chiral partner, the sigma. One of the 
goals of this work is to develop practical methods to extract critical indices from chiral symmetry breaking 
transitions and classify critical points as Gaussian (trivial) or non-Gaussian (nontrivial). We shall see that 
the ratio R„ = /M^ and its dependence on coupling and bare fermion mass predicts the critical coupling 
and the the critical index 5. We shall relate i?,r to the transverse and longitudinal susceptibilities of the chiral 
phase transition and show that R.^ is completely determined by the theory's equation of state. The resulting 
formula has already proved its usefulness in studies of lattice QED [N=2] and we anticipate applications 
to the other models listed above. Most past work on lattice QED and its chiral symmetry breaking phase 
transition have concentrated on the order parameter < ipijj > and its equation of state [7]. Since these 
quantities are measured at nonvanishing bare fermion mass, the critical point must always be inferred and a 
rapid crossover between a symmetric and a broken phase occurs in raw simulation data. This obscures the 
critical point and the value of critical indices. However, since is determined by the same equation of state, 
any hypothesis one makes concerning < > and its equation of state must predict or be discarded. 
This fact should help us find the correct theoretical ideas describing the phase transition of lattice QED, 
and should be applicable to other models as well. 

Following the same philosophy, we shall study as a function of < -tjjil) > itself. We shall see that 
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combining correlation length scaling with the Goldstone description of the pion, yields particularly effective 

methods to find the critical coupling and the anomalous dimension 77. Since and < '4)xp > are the easiest 
observables to measure in simulations and typically the least technical quantities to calculate in models, this 
approach to rj shoiild prove quite useful. In addition, the anomalous dimension rj is central to discussions of 
nontriviality and thus lies at the heart of the issue of whether non- asymptotically free theories exist. 

The rest of this paper is organized as follows. In Sec. 2 we introduce our general theoretical framework and 
discuss renormalization group trajectories, the chiral equation of state and the physics behind hyperscaling. 
In Sec.3 we concentrate on ii^ and show that it is determined by the chiral equation of state and can be 
used, because of the Goldstone nature of the pion, to locate the critical coupling and measure the critical 
index 5 in an (almost) model independent fashion. Sec. 4 considers the chiral equation of state and the 
constraints that PCAC places on it. In Sec. 5 the relation of the pion mass, and the chiral condensate is 
developed resulting in a novel way to find the critical point and the anomalous dimension r]. Sec.6 considers 
the fermion mass and the index S and discusses the scaling region where the analysis of this paper applies. 
It also discusses non-asymptotically free dynamics and the breaking of scale invariance at the critical point. 
In Sec. 7 we illustrate some of our analyses in the context of two models, the four-fermi and the linear sigma 
model. Sec.8 contains some brief concluding remarks. 

2. RENORMALIZATION GROUP TRAJECTORIES AND HYPERSCALING 

The general idea of renormalizability is that a theory's cutoff dependence can be absorbed into a finite 
set of bare parameters in such a way that its low energy physics is insensitive to the cutoff. Once this is 
done, it is possible to find the lines of constant physics (RG trajectories) in the bare parameter space. These 
lines are uniquely defined no matter what observable is taken. In this way low-energy quantities depend on 
each other and not on the cutoff. We will see that this is possible if hyperscaling is obeyed. 

In our discussion we shall assume that the reader is familiar with traditional statistical mechanics topics 
of homogeneity, the Equation of State (EOS) of ferromagnets, correlation length scaling etc. We shall pass 
between the languages of field theory and statistical mechanics freely in much the same spirit as the textbook 
by C. Itzykson and J.-M. Drouffe [8]. 

The hyperscaling hypothesis claims that the only relevant scale in the critical region is the macroscopic 
correlation length ^. If this hypothesis holds, it is possible to do dimensional analysis using this correlation 
length as a scale. As a consequence, all dimensionless observables, e.g. mass ratios and renormalized 
couplings, should be independent of The hyperscaling hypothesis can be stated as 

Fsing = t'-'^Fih/t^) ^ (2.1) 

where notation is standard (F is the free energy, h is the external magnetic field, t is the deviation of 
the dimensionless temperature (coupling) from the critical point, etc.). All thermodynamic quantities are 
obtained by taking the derivatives of the free energy. In particular, the order parameter, defined as < ^ >= 
dFsing/dh, satisfies the Equation of State (EOS) 
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<(j)>= t^-''-^F'{h/t^) = t^F'{h/t^) (2.2) 

This relation also defines the magnetic exponent: (3 = 2 — a — A. Similarly, the susceptibility exponent, 7, 
is obtained from 

^ = = t^-^F"{h/t^) = t-'^F"{h/t^) (2.3) 

i.e. 7 = A - 0. 

For a given channel, the corresponding correlation length (= inverse mass) is defined by 

_ L 1^1' < 'Pi^)m > 
'^^^ - L < 4>{x)m > ^'-'^ 

If the field develops a vacuum expectation value, then connected correlation functions should be taken in 
eq.(2.4). This will be understood throughout the paper. In the scaling region its behavior is given by 

i = t-''g{h/t^) (2.5) 

Combined with eq.(2.1), it leads to the hyperscaling relation between the critical exponents: di' = 2 — a. The 
other relations between the exponents follow from the scaling form of the two-point function: < (j){x)(j){Q) >= 
l/|xp'^*/(x/C) [9]. 

Since hyperscaling is an important concept, we should explain its meaning and outline possible impli- 
cations. It is generally believed that violations of hyperscaling lead to triviality. This is due to various 
inequalities between certain combinations of critical indices [10]. In simple models like scalar field theories 
and spin systems, the quantity that measures the violation of hyperscaling is the dimensionless renormalized 
coupling [11]. It is defined through the nonlinear susceptibility x^"'^ 

^« = -^. X'"> = ^f- (2.6) 

The renormalized coupling is essentially the properly normalized connected four-point function. In terms of 
correlation functions, the nonlinear susceptibility is the zero-momentum projection of the four-point function 



("0 = f < Ct>{x)4>{y)ct>{z)m >conn (2.7) 
J xyz 

The term in the denominator of eq.(2.6) is to account for an extra integration. By differentiating the free 
energy and using the EOS for the correlation length (eq.(2.5)), we arrive at 

gji = t-^^+^+'^''H{h/t^) (2.8) 
We can trade t for the correlation length in eq.(2.4) and rewrite eq.(2.8) in terms of ^ as 



= ^(2^-T-'''')/-^(/l/f^) 
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(2.9) 



In specific (ferromagnetic) models an estimate of the exponent can be made: since multi-spin correlations 
can not extend over a larger range than pair correlations, the following inequality holds [10] 

2A<7-|-(iz/ (2.10) 

Because qr is dimensionless, hyperscaling implies that it is independent of ^ and the exponent in cq.(2.9) 
must vanish. In this case the renormalizcd coupling is a function of only one bare variable. Clearly, strict 
inequality in eq.(2.10) implies triviality. We mention that an equivalent way of stating the above inequality 
IS dv >2 — a [12]. It means that the singular part of the free energy vanishes at most as fast as required by 
hyperscaling i.e. Fging > 

In a similar fashion the scaling of the mass ratios can be derived [13]. If hyperscaling is satisfied, the 
ratio of any particular pair of masses satisfies 

R{t,h) = G{h/t^) (2.11) 

where G{y) is a universal function. Comparing it with the expression for the renormalizcd coupling, we see 
that both observables depend on the same variable. One of the relations can be inverted to solve for the 
bare variable e.g. h/t^ = H~^{gji). This defines an RG trajectory for each value of qr and can be used to 
obtain the relation between the two observables R ~ R{gR,)- Note that this relation is independent of the 
bare parameters. The same manipulation can be done with two mass ratios. We note that the important 
point in the inversion is that both observables depend on just one bare variable, so that the inverse relation 
can always be found, at least in some regions of parameter space. 



3. a-TT SPLITTING, CRITICAL EXPONENTS AND EQUATION OF STATE 

In this section we will relate the a — tt splitting with other universal quantities like critical indices and 
universal amplitude ratios. Since the scalar and pseudoscalar propagators represent fluctuations of the order 
parameter, information about the ct — tt spectrum is contained in the equation of state. Before we consider 
the derivation, we briefly illustrate the connection between the behavior of < •ip'ip > and the a — tt splitting. 
We use the spectral representation to express the order parameter 

<^i»=j_^ .Ap(A)^^ (3.1) 
where p{X) > is the eigenvalue density of the Dirac operator. We define the two susceptibilities 

Xa= < ^V(a;)V'V'(o) >, Xtt = / < iph5i'ix)il>h5i^{0) > (3.2) 

which are the zero-momentum projections of the scalar and pseudoscalar propagators. Again, we emphasize 
that the scalar correlation function is understood as the connected one. These two susceptibilities are 
related to the order parameter via: Xa = d < -tjjil) > /dm and the Ward identity Xn =< V'V' > /m. (We 
will derive this relationship later in this section.) With this in mind, it is straightforward to show that these 
susceptibilities can be rewritten in terms of the spectral function as 
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X. = X.-y_ dXpiX)^^^-^^ (3.3) 

The integral on the right hand side is positive leading to the following inequality Xct < Since the mass, 
defined as in eq.(2.4), is related to the susceptibility via = Zx~^, and since scalars and pseudoscalars 
renormalize in the same way, we obtain the following inequality (level ordering): < M^. 

From cqs.(3.1) and (3.3) it is apparent that both, the order parameter and the (t — tt splitting disappear 
in the chiral limit if the zero-mode is absent from the spectrum i.e. when p(0) = 0. Conversely, the presence 
of the zero-mode induces both a nonzero order parameter, < ip^p >— TTp{0) ^ 0, and cr — tt splitting. The 
same mechanism that leads to a nonvanishing value of the order parameter is responsible for lifting the 
degeneracy within the parity doublets. 

The starting point in our analysis is the scaling form of the equation of state [8,9]. In the case of 
continuous symmetry the familiar form h = f[t/M^/^) becomes 

K = M^M^-^f{t/M^I^) (3.4) 
where, for simplicity, we denote the vector order parameter and external field by Mq and ha respectively. 



and the modulus of the order parameter as M = ^/M^M^. Our convention is that t > corresponds to the 
broken phase. Eq. (3.4) is the gap equation for the order parameter. In the limit when h ^ we have 

= M^f{x) (3.5) 

which has a nontrivial solution M ^ only if the function f{x) has a positive zero f{xo) = 0. The 
spontaneous magnetization is given by M — {t/xo)^- This determines only the modulus of the order 
parameter (vacuum degeneracy), whereas the orientation and, thus, the particular Hilbcrt space is fixed by 
the direction of the external field. The function f{x/xo) = f{x)/f{0) is the same for any theory within a 
given universality class. At the critical point a; = the response of the system is singular and is given by 
h = M^f{0) > 0. In order to have no spontaneous magnetization in the symmetric phase f{x) should be 
positive on the negative x-axis. 

The response to an external field ha is given by the inverse susceptibility {x~^)ab = dha/dMh which 
can be obtained from the equation of state 

ix-'U = Sa,M^-'m + {5- l)^M^-V(x) - ^^M^-V'(x) (3.6) 
where x = t/M^/'^. Eq.(3.6) can be rearranged into 

(X-^. = (Sa, - ^) M^-^fix) + ^M-^ (Sfix) - ^/'(x)) (3.7) 
in order to separate the susceptibility into transverse and longitudinal parts 



Xt' = M'-'f{x), xl' = M'-' (^Sfix) - ^f{x)^ 



(3.8) 



The expression for the transverse susceptibility is the Ward identity which can be rewritten, after using the 
EOS, as 
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h = Mxt' (3.9) 

In the broken phase where M ^ 0, the divergence of the transverse susceptibihty in the /i — > Hmit signals 
the appearance of massless modes (Goldstone bosons). At the critical point {x = 0), eqs.(3.8) imply that 

= (5 (3.10) 

Xt 

This defines the exponent (5 as a measure of the relative strength of longitudinal and transverse responses of 
the system at the critical point. Eq.(3.10) can also be derived by noting that, at the critical point, M ~ /i^/*, 
and using xl = dM/dh and xt = M/h. 

When the symmetry in question is chiral symmetry, eqs.(3.8) can be combined into a statement about 
the meson masses. The external field is replaced by the bare fcrmion mass m and the order parameter, 
Ma, has components (< ip-ip >, < iJji"fi:,T'^'ip >) for the flavor group whose generators are T". ( The correct 
mapping to ensure dimensional consistency is < ^ip >—> K^M and h —> mK^ . In the following we will work 
with dimensionless quantities, so we are allowed to take A = 1 from here on). 

We introduce the bare mass in the standard way through rmpip, so that the ground state is parity 
invariant. Then, condensation occurs in the scalar channel and the only nonvanishing component of the 
order parameter is < -tjjtp >. The EOS, eq. (3.4), then reads 

TO =< V^V fit/ < V'V' >^^'^) (3.11) 

The sigma and the pion are longitudinal and transverse modes, respectively. The Ward identity, eq.(3.9), 
now becomes 

< >= TO / < ■il)i^5il){x)il)i^5ip{0) > (3.12) 

The masses are related to the corresponding susceptibilities via — Z^^x^^, = Z^Xa^- Because of 
chiral symmetry, the wavefunction renormalization constants for the two modes are equal and the mass ratio 
(squared) is the susceptibility ratio i.e. 



Thus, eqs.(3.8) implies 



= Ml/Ml = Xt^Ix'l (3.13) 



^ =5_|Z^, ^^(0,m) = ^ (3.14) 



i?^(t,m) /3 !{x) ' ' 5 

Eq.(3.14) is a stronger statement than the Goldstone theorem eqs.(3.9) or (3.12). Not only does it contain 
the above Ward identity, but it determines the splitting within the parity multiplet away from the chiral 
limit in both phases. Before we discuss the above result, we recall the EOS for the masses, 

M„{t,m) = i^5a(^), M^{t,m) = f^g^j^) (3.15) 

While the equation for M„ follows from hyperscaling, the same reasoning does not lead to the analogous 
expression for Its form is fixed, however, by eq.(3.13) which states that the ratio R„ is a function of x and. 
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thus, of m/t^, implying eq.(3.15). The functions g{y) are universal up to a multiplicative constant. Clearly, 
the Goldstone nature of the pion imposes a different boundary condition on requiring ^^(y) ~ y^^"^ 

near the origin. The argument m/t^ is understood in the sense sgn{t)\t\^ so that it changes sign 

as one goes through the critical point. This insures single- valucdncss of g{y) and the possibility of global 
inversion. Taking the mass ratio eliminates non-universal multiplicative factors as well as the exponent v 
making it a universal function oi y = m/t^ only 

ii.(t,m) = ^=G(^) (3.16) 

In this equation, we have the low energy observable on one side and the bare parameters on the other. From 
eq.(3.14) it follows that 

R^{t = 0, m) = G(oo) = I (3.17) 

There are two limiting values of i?,r that are fixed by chiral symmetry. In the broken phase {t > 0) pions 
are Goldstone bosons, so i?^(f > 0,m = 0) = 0. In the symmetric phase the sigma and pion degeneracy 
implies R^{t < 0, m = 0) = 1. The ratio can never exceed unity since Euclidean propagators satisfy 
|i^7r(a;)| > |r>cr(a;)| which, in the large |a;| limit gives < 1. As the bare mass increases, the ratio approaches 
the same value regardless of which phase it originates from. This is so simply because, as one increases the 
amount of explicit symmetry breaking, the dynamics becomes insensitive to the type of symmetry realization 
in the vacuum. Thus, if we fix t and plot the curve R^rit, m) versus m (Fig.l), in the broken phase i?^ = i?> 
will increase from zero while in the symmetric phase R^^ — i?< will decrease with m starting from i?7r = 1. 
Both families of curves approach R^ = l/S asymptotically from above (symmetric) and below (broken). 
Because of the scaling form of eq.(3.16), Rj^it, m) = G{m/t^), it is clear that Rc = G{oo) is independent of 
m. So, we have traded small t for large m. Since < R„ < 1, it follows that R> and R^ curves have slopes 
with opposite signs. Thus, all the curves in the broken phase are monotonically increasing and lie below 
R = Rc i.e. Ry{m) < Rc and R'^{m) > for all m. Similarly, in the symmetric phase R<{m) > Rc and 
i?< (m) < 0. Therefore, the following inequality holds 

R>{m) < Rc < R<{m) (3.18) 

which is saturated asymptotically for large values of m. Since Rc = l/S, eq.(3.18) gives a bound on the value 
of (5 

^ < ' < ^ ^'-''^ 
Any curve Rj^lrn) in the broken phase produces an upper boimd on S etc.. This bound improves as the 
ratio m/t^ becomes larger. The analysis of the mass ratio i?^ is consistent with the function /(.x) being 
positive semi-definite with only one zero at xq (Fig. 2). A possible change in its monotonicity would imply 
the vanishing of its first derivative for at least one non-zero value of x. This would be in conflict with the 
physical behavior of R„. 

There are several reasons why this result is useful. Its application to data analysis produced by lattice 
simulations is obvious. Recent applications of this method to strongly coupled QED have been very successful 
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[14], giving more accurate estimates of the critical coupling and the exponent 6 than those obtained with 

conventional methods on data samples of comparable statistics. 

The reason for this can be easily understood. Instead of dealing with extrapolated data, as normally 
done when studying chiral symmetry breaking through < ijjip >, here we determine both, the critical coupling 
and critical exponent S from the raw data. Furthermore, the scaling form for the mass ratio is more accurate 
for larger values of m/t^. This means that, instead of simulating at the critical point and looking at the 
m — > limit, we can work away from it and use large masses without losing accuracy. The key point, from 
the point of view of numerical simulations, is the m-independence of Rj^it, m) at the critical point, eq (3.17). 
This can exploited by plotting R„{t,m) as a function of f3 for different m values. The spectral method can 
be used to get very accurate plots, which cross at the critical point, thus determining /3c. This is reminiscent 
of the techniques of Finite Size Scaling when applied to the Binder cumulant, as noticed also by Boyd et al. 
in ref. [15] where a similar method has been proposed and successfully applied to the QCD chiral transition 
in the strongly coupled phase. 



4. PCAC, HYPERSCALING AND TRIVIALITY 

Eq.(3.14) is a differential equation that determines the function f{x) in terms of other universal quanti- 
ties like critical exponents and mass ratios. Its solution requires one boundary condition which can be fixed 
at an arbitrary value of x. The formal solution for f{x) is given by the form 



/(x) = /(0)exp(/3 r^{6-±) 
\ Jo y -n-TT 



(4.1) 



which shows that the ratio f{x)/f{0) is universal. In the symmetric phase, a; < 0, — > 1 in the chiral limit 
{x — > — oo). Therefore, in this limit eq.(3.14) gives 

1 = P{6-1)= hm (4.2) 

X— >-(X) J[X) 

where we used the scaling relation between the critical exponents: 7 = f3{S — 1). The above equation implies 
that, for X — > —00, f{x) w C{—x)'^. In general, the expansion of f{x) in the symmetric phase has the form 

Another constraint on the behavior of f{x) is PCAC. It refers to the chiral limit in the broken phase. 
Eq.(3.14) can be rewritten as 

Pfjx) 

'''^ S0fix)-xf'ix) ^^■"'^ 

which implies the vanishing of the pion mass in the chiral limit (a; — > xq). Now, we show how PCAC forces 
xo to be a first order zero. Assume that around xq, f{x) vanishes as 



fix) « a{xQ - xY 
In the chiral limit (in the broken phase) eq.(4.3) becomes 
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(4.4) 



{6l3- p){xo - x) + pxo 
so vanishes linearly. For x xo the equation of state is 

= fix) « a{xo - xY (4.6) 

Clearly, PCAC, which requires m, forces p = 1. Substituting eq.(4.6) into the expression for 

(eq.(4.3)) gives 

^ _ Q ^^1^^^ u 7) 

{5(3 - l)m + axo<il)tp>^ ^ ' ' 

As m — > the leading contribution is 

axo <-ipip > 

When compared with the PCAC relation, = 2m < 'iptj; > / f^, eq.(4.8) gives 



^^-£-^TrfT7 (4-8) 



M2 p M4 ^ ' 

This relates the pion decay constant to the physical scale. We can use the definition of the critical exponents 
to study the behavior of in the scaling region e.g. P{S + 1) = 2A — 7. However, using the EOS for < ?/)^ > 
(eq.(3.11)) and the EOS for M„ (eq.(3.15)) the above relation can be written as 

a 

with K{y) being a universal function. The exponent is the same one that appears in the expression for the 
renormalized coupling (eq.(2.9)). Clearly, since the above ratio is dimensionless, hyperscaling implies the 
vanishing of the exponent. Conversely, its violation would imply that the above ratio is cutoff dependent. In 
theories with Yukawa couplings, triviality is intimately related to the issue of compositeness of the scalars. 
For a nontrivial continuum limit to exist, it is necessary for fermions to exchange composite scalars. This 
is contained in the Goldberger-Trciman relation g = Mp/f-^, which, after recognizing that the pion radius 
scales as r^r ~ I/Ztt implies that the coupling vanishes in the limit where pions are pointlike (g ^ Mpr.^). 
In that sense, compositeness means that the pion decay constant has to scale as the physical mass scale. 
Pointlike structure, on the other hand, implies that diverges in physical units. Although Baker's inequality 
is rigorously proven only for ferromagnetic systems [10] , it is interesting to note that it is valid in the above 
equation in all known models where f-^jM^ either diverges or approaches a constant depending on whether 
pions are pointlike or not. The inequality A < 7 + 4z/ is in agreement with it. 



5. PION MASS, ORDER PARAMETER AND ANOMALOUS DIMENSIONS 

Nonvanishing anomalous dimensions and compositeness are intimately related. In addition, the compos- 
iteness condition itself is tied to the existence of a fixed point of the underlying theory. It is reflected in the 
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vanishing of the wave function renormahzation constant associated with the composite degrees of freedom. 
Let us explain this in some detail. Consider a general field theory where a bound state \B >, with binding 
energy Eb < 0, appears. From general considerations [17] , the wave function renormahzation constant Z 
can be expressed as 

Z = J2\<b,E\B>\^ (5.1) 

b,E 

where \b,E > is the bare (elementary particle) state with energy E. Using standard techniques from quantum 
mechanics, it can be shown that Z satisfies the following equation [17] 

with C^iE) > being the total decay rate of the state \B >, C^iE) = 2'Kj2b I < b,E\B > j^. As such G{E) 
is proportional to an effective coupling constant. The compositeness condition, Z = 0, implies two things: 
1) The composite state has no projection in the space of bare states (eq.(5.1)) i.e. < b\B >= for any \b >, 
and 2) It is a sum rule that places an upper bound on the effective coupling, G^{E), (eq.(5.2)) and can be 
interpreted as a fixed point condition. 

As the fixed point is approached, the order parameter, wave function renormahzation constant, Z, 
and all the masses vanish. If there is a single correlation length, the vanishing of Z is determined by the 
anomalous dimension, r], via Z ^ Also, the order parameter, < (f) > scales as < (p >'~ C"***) where 
= {d — 2 + r])/2 is the scaling dimension of the field (j). Combined with the standard scaling laws in the 
symmetric limit, < (p >~ and (, t~'^, this leads to the following scaling relation between the exponents 

^ = l{d-2 + v) (5.3) 

The appearance of the anomalous dimensions guarantees the compositeness of the degrees of freedom 
involved, a necessary condition to produce a nontrivial continuum limit. Consider the pair of equations used 
in previous sections 

m =< >^ f{x), ' =< V^V' f{x) (5.4) 

We shall sec how eqs.(5.4) restricts the dependence of on the order parameter in the chiral limit 
and show that this behavior is universal and is determined by the magnitude of the anomalous dimension 
rj. As such, it contains information about the continuum limit of the theory and is capable of distinguishing 
between mean field and non-mean field behavior. This being the case, perhaps, the best way to start is to 
recall this dependence in the cr-model in four dimensions. The masses and EOS are given by 

Ml = -t+^v\ Ml = -t+\v\ h = v{-t+^v^) (5.5) 

where v =< cr > is the order parameter. The universal function in this case is a straight line f{x) = —x+X/6. 
In the chiral limit the order parameter is obtained from f{xo) = giving Vq = 6t/X. Since this is a mean 
field model, there is no wavefunction renormahzation and = x~^. Eq.(5.5) thus, reproduces the Ward 
Identity and EOS of eq.(5.4). For a fixed value of t, is a linear function of v"^ with different intersections 
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depending on the phase . The expression for the pion mass can be rewritten as = (A/6)(w^ ~ ^o)- Thus, 
as discussed in the previous section, the pion mass vanishes hnearly as the chiral hmit is approached . In the 
symmetric phase the v hmit resuhs in a nonvanishing pion mass . Again, the chiral hmit is approaclicd 
hnearly. The reason for such a simple behavior is the absence of wavcfunction rcnormalization (vanishing 
anomalous dimensions). It is easily demonstrated that, as long as the identification = can be made, 
eqs.(5.4) insures the linear dependence between and < -tptjj To show this, we recall the expansion 
of the universal function f{x) around the chiral limit in the symmetric phase: f{x) ~ which leads, in 
the limit < -tjjil) >— > 0, to \t\'^ . The general dependence of the pion mass on the order parameter for a 

mean field theory is summarized in Fig. 3. 

Since, for nonvanishing anomalous dimension r? > 0, the wavefunction renormalization constant scales 
as well, it will give some curvature to near the origin both in the symmetric phase and at the critical 
point (Fig. 4). This curvature will be a signal of nontrivial behavior of the theory and is related to the 
compositencss of the pious. This curvature is 1) determined by the magnitude of the anomalous dimensions 
and 2) is opposite from that induced by finite size effects. In this way, the above plot contains information 
about both the thermodynamic limit and nontriviality of the theory. 

In general, the pion mass is related to the susceptibility as = Zx~^. The wavefunction renor- 
malization constant scales as Z <^ and, since the order parameter scales as < -ipil) >'~ we have 
Z ~< VSV) >'"'/'^. At the critical point, t = 0, ~< '^V' >''"^- Thus, 

~< >^-i+W/5 (5.6) 

The exponent can be expressed in a slightly difl[erent form using the scaling relations between the critical 
exponents: (3{6 — 1) = 7 and 7 = u{2 — r]). This leads to 

In four dimensions, the expression in parenthesis is 1/(1 + ??/2) < 1. Thus, at the critical point we have 

Ml - (< >'^Y"^, {t = 0) (5.8) 

The important point is that the slope of the curve is infinite near the origin because of the anomalous 
dimensions: (S/v = l + r//2 in four dimensions. 

In the symmetric phase, we use the asymptotic expansion of the universal function f{x) for x —00 

[16]: 

00 

fix) = ankp-'"^ « aolxp + a,\xr^^ (5.9) 
The pion susceptibility and its mass are determined from this expression giving 

m; ~< ^ («,^J^1^ + «7dSS«7?) 

The first term gives the leading contribution near the origin. To obtain its magnitude, we use once again 
(3(6 — 1) = 7. In that case the leading term gives 
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Ml ~ (< >2)W2/3^ < 0) (5.11) 
Prom the scaling relations it follows that the exponent in eq.(5.11) is less than 1: 

!^ = 1- <1 (5.12) 

2(3 (1-2 + 1) ^ ' 

Clearly, for vanishing r], the leading term is constant and the second term in eq.(5.10) determines the 
curvature. The exponents in that case are given by mean field theory, so 1//3 = 2. Therefore, approaches 
some finite value linearly. The curvature, and the infinite slope come from the anomalous dimensions. 

Finally, in the broken phase, t > 0, the order parameter is nonvanishing in the m — > limit and the 
function f{x) vanishes for x = xo. PCAC constrains this to be a first order zero i.e. f{x) « a{xo — x). 
Therefore, in the chiral limit we have 

Ml = Zx~^ ~< lAV' >^-^+'^'''^ a{xo - x) ~ (< V^V > - < V'V' >o), {t > 0) (5.13) 

and the pion mass vanishes linearly. 

A remark about the apparent vanishing of the pion mass in the symmetric phase (in the chiral limit) 
should be made. At first glance this result is surprising, because symmetry considerations alone do not force 
its vanishing. In mean field theory where — everything is canonical and equation (5.10 ) tells us 
that Ml is finite in the symmetric phase even for m = 0: 



Xtt 



1 ~< >*-^ {\t\/ < vv- >^^'^r ~ i*r (5-14) 



The vanishing comes from the anomalous dimensions i.e. wavefunction renormalization! Let us analyze this 
point in more detail. The correlation length in Euclidean theory is 



.2 _ /. \^\'D{x) 



or, in momentum space 



This is a standard way of extracting the infra-red piece of the propagator that dominates at large separations. 
If we denote the self-energy corrections to the meson propagator as n(fc^), the full propagator in momentum 
space reads 



£)-l(fc2) = ^2^ ^2 ^ jj(^2) ^ ^2 ^ ^2 ^ jj(Q) ^ fc2n'(0) + Ilik'^) (5.16) 

This decomposition is made so that U{k'^) = 0{k'^) and the first two terms dominate the infra-red region. 
After identifying the fc^ term in eq.(5.15b) as the inverse wavefunction renormalization constant, = 
1 + n'(0) and inverse susceptibility as = Mq + 11(0), the meson propagator can be written as 

D-\e) = Z-\e + Zx~^ + nij(A;2)) « Z-\k'^ + M^ + 0{k^)) (5.17) 
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Thus, at low momenta, the propagator resembles that of a meson with mass = Zx~^- This is the 
definition of the mass as given by eqs. (5.15). In this case analytic continuation to Minkowski space is 
standard and the pole structure is recovered. As is obvious from the decomposition in eq.(5.16), such 
manipulations assume a particular analytic structure of the composite propagators which is guaranteed only 
if the fermions are massive. Additional nonanalyticities in the form of branch cuts appear when fermions 
are massless as occurs in the chiral limit in the symmetric phase. The anomalous scaling of the scalar 
propagator is a consequence of these nonanalyticities - the fc^ term is absent. Rather, the leading low- 
momentum behavior is given by 

The derivative diverges for fc^ = giving zero mass as defined by eqs. (5. 15), although there might be a pole 
in Minkowski space. Such a propagator does not have an exponential, but a power law decay. So, there are 
long-range correlations ("massless modes"), but the theory is not scale invariant. (Clearly, the case r] = 
is well behaved. In the broken phase there are no infra-red "problems" of this sort because fermions are 
massive and eq.(5.15) is the appropriate definition of the correlation length.) Some attempts to confront 
this problem in a lattice study of three-dimensional Gross-Neveu model have been reported in [4] . 

When studying the theory in a finite volume, as is always done in lattice simulations, one has to insure 
that results are not obscured by finite size effects. Since different quantities have different sensitivity to 
finite size effects, the effective thermodynamic limit, therefore, must be monitored carefully through the 
consistency of certain relations. One set of such relations consists of the Ward Identity and the Equation of 
State. Being the lightest particle in the spectrum, the pion is most sensitive to finite size effects. On a finite 
lattice, the pion mass would tend to a higher value than in the thermodynamic limit. This would result in 
a change of sign in the curvature. In that sense, the plot vs < ■ip'ip >^ is also suited for controlling finite 
size effects in theories with nontrivial fixed points. 



6. HEAVY QUARK LIMIT AND EXPONENT d 

We have seen that = G{m/t^) with G{oo) = 1/6. This means that, in the large-m limit, R„ — > 1/6 
(not necessarily 1). At first sight this looks puzzling since one expects that the nonrelativistic limit can 
be taken for very large constituent masses so that any meson mass should approach twice the constituent 
mass i.e. Mmeson ~ 2m which would necessarily imply i?^ k, 1. We should clarify what is meant by a large 
fermion mass: the bare mass is always much smaller than the cutoff and large m means large compared to a 
typical mass scale, for example M^r in the chiral limit. In this way we are sure to stay in the scaling region 
where universality arguments hold and the RG trajectories are \miquely defined. 

In theories that are not asymptotically free expansions in powers of p/m are poorly behaved because 
the force between constituent quarks is strong at short distances forcing a large kinetic energy due to 
uncertainty relations. Thus, whatever the quark mass, it is never large compared to a typical kinetic energy. 
The composites in the scaling region are always relativistic. We note also that the bound states in the 
two phases are quite different from each other. In the symmetric phase, an increase in the bare mass is 
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compensated by a decrease in the coupling in order to keep the ratio unchanged. This is expected because 
the zero-point energy is reduced by increasing the mass and less attraction is needed to produce the same 
effect. In the broken phase, the opposite happens. 

In theories hke QCD this is not so because the force between tlie quarks is weak at short distances and 
the typical momenta could be small compared to the bare mass. Thus, in principle, one can still be in the 
scaling region and have a heavy quark limit. The most direct way of understanding the differences between 
asymptotically free and non-asymptotically free theories is to note that in the scaling region of the former 
\t\ <C 1 means weak coupling {qc = 0,t = g^), whereas for the latter \t\ <C 1 means strong coupling. 

We should add that ^ = 1 has different physical origins in asymptotically free and non-asymptotically free 
theories. Since the exponent 5 gives the response of the system at the critical coupling: m =< V'V' \t=o-, 
in asymptotically free theories the system responds as a free theory, so 5 = 1. In that context, Rt^ « 1 is 
a reflection of the fact that there is very weak binding which could not possibly distinguish between the 
scalars and pseudoscalars. In strongly coupled non-asymptotically free theories, tightly bound relativistic 
composites are formed in the scaling region. Their presence at high energies is the main difference relative to 
QCD-like theories where the only relevant degrees of freedom in the scaling region are quarks and gluons and 
where binding occurs in the infra-red regime. Because of the non-asymptotically free nature of the couplings, 
the ultra-violet asymptotics of the scalar correlation functions is not canonical - the theory has anomalous 
dimensions. In terms of the anomalous dimensions the exponent 5 in four dimensions is 

5 = 1^ (6.1) 

Here, ^ = 1 is a consequence of the large anomalous dimension, r] = 2. The meaning of this particular limit 
can be best understood if we write down the first two terms in the operator product expansion of the fermion 
propagator [18] 

« ^ + ^ + + - (6-2) 

As 77 ^ 2 i.e. ^ ^ 1, the system reacts to the bare mass the same way it would react to a change in the 
dynamical mass of its constituents. Thus, the persistence of R„ = 1 away from the chiral limit in this case 
means that the system can not distinguish between bare and dynamical masses. 

As a final remark, we use once again eq.(3.14) to argue the absence of the dilaton and scale invariance 
in theories with spontaneously broken chiral symmetry. Unlike previous treatments [19] our argument will 
require no knowledge of the composite propagators and is thus independent of the approximation scheme. 
Instead it follows from the scaling form of the EOS only. The idea of a dilaton as a Goldstone boson of 
spontaneously broken scale symmetry was introduced some time ago in ref. [20] . If the theory is scale invariant 
and if it breaks chiral symmetry spontaneously, it, at the same time, generates a scale, the fermion mass, 
and is no longer scale invariant. Since scale invariance is thus broken "spontaneously", one naively expects 
that there should be a corresponding Goldstone boson, a massless scalar, in the spectrum. 

We use the scaling relation /? = 7/(5 — 1) to rewrite eq.(3.14) as 

= 5-x{5-l)/^{f'{x)/f{x)) ^^-^^ 
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In general, it is clear that there has to be splitting between the a and tt i.e. there can be no massless scalar 
in the theory as long as chiral symmetry is realized in the Nambu Goldstone manner. An exception might 
only be the theory with 5 = 1. In that case R^^ = 1 for any value of the bare parameters - not only are 
the parity partners degenerate, but they respond to symmetry breaking in the same way. This would be a 
very unusual realization of chiral symmetry: a and tt are indistinguishable for any finite m. In QCD, where 
6=1, this situation is avoided by explicit violations of scaling due to quantum corrections. These scaling 
violations are the sole source of interaction and of the a — it splitting. In other cases e.g. strongly coupled 
QED with 6 = 1, the a — w degeneracy in the broken phase could not be reconciled with PC AC which 
requires i?^ m. Thus, if the current algebra relations are to be realized, scale invariance must be violated. 
The entire content of PCAC is contained in these scaling violations. Consequently, the mass of the a comes 
solely from the scaling violation. 

7. TWO EXAMPLES (SCALING PLOT) 

To illustrate and complement the previous discussion, we analyse two simple examples: the four-fermi 
and linear a models. In addition to demonstrating realizations of the general ideas in these two models, we 
also discuss the mutual dependence of the mass ratios , or scaling plots, (as introduced in lattice spectroscopy 
in ref. [21]), and the heavy quark limit in these models. The scaling plot will prove especially suitable to 
argue the equivalence of the two models. (For notations and details of computations related to this Section 
see [13].) 

First, we start with the four-fermi model. In the leading order in 1/N there is no distinction (on the 
technical level) between discrete and continuous chiral symmetry. We consider the continuous U{1) x U{1) 
model. For 2 < d < 4 these models are renormalizable and nontrivial. To leading order, the critical exponents 
are f3 = i' = l/{d — 2),6 = d— l,r] = 4: — d. The gap equation for the fermion self-energy S is given by the 
tadpole contribution 

T, = m + 4:g^ <i)^> (7.1) 

or, explicitly, in terms of S, 

c5^E2-=t+| (7.2) 

where e = 4-d (not necessarily small, < e < 2) and c = 46/(2 - e), with b = 2r(e/2)/(47r)''/2_ -y^g 
express all the quantities in units of the momentum cutoff A = 1. The masses of the composites M„, are 
given by the one-loop diagrams [13] 

= ^Z, M^ = M^+ 4S2 (7.3) 

with the wave function renormalization constant Z = YJ^ jhg^ . 

We first discuss the dependence of on the order parameter. In this case, the expectation value of 
the scalar field can be used to make contact with the cr- model i.e. S = g < tr >. Clearly, the relations in 
eq.(7.3) can be combined to give 
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Ml = ^i:\cg^i:''-^-t) (7.4) 

which, in the Hmit = (rf = 4) reduces to the cr-model expression (eq.(5.5)). In particular, at the critical 
point, t = Q, S^. After recognizing that (3 = u in this model, we recover eq.(5.8). Similarly, in the 

symmetric phase, near the origin and for fixed i < 0, we obtain ~ which is eq.(5.11). Thus, all 

the general features of this plot arc clear from eq.(7.4): linear behavior in the broken phase, concavity in the 
symmetric phase and the vanishing of the pion mass in both phases. The last point is especially interesting 
considering the discussion of this problem in sect. 5. One can illustrate this point further by calculating the 
pion propagator in the chiral limit in symmetric phase. Its form is given by [4] 

D-\k^) = -t + Ck''-'' (7.5) 

which has a pole in the complex fc^-plane (in Minkowski space), but gives a power law behavior of the 
Euclidean correlator. 

As we remarked before, the general idea behind the universal behavior of a particular dimensionless 
observable is that, given its functional form R = G{m/t^), one can invert this relation to find an RG 
trajectory and express another observable as R' = R'{R). We define two mass ratios: R^ = M^/M^ and 
Rp = 4Mp/M^ in terms of which the above equation becomes 

R^ = l- Rf (7.6) 

In the chiral limit R.^ = and Rp = 1 i.e. = 2Mp in the broken phase, whereas in the symmetric phase 
i?7r 1 and Rp ^ as m ^ 0. The curve i?^ versus Rp is universal. The prediction of the four fermi 
theory (eq.7.6) is a straight line with unit slope. All the points on the curve that lie below i?^ = 1/6 belong 
to the broken phase and those above to the symmetric phase. All the physical points lie on the straight line 
eq.(7.6) . Clearly, the naive heavy quark limit would require R^ = Rp = 1 which is completely missing from 
the plot. The points below the line R^ = 1/6 are in the broken phase and at the critical point Rp = 1 — 1/6. 
These two values correspond to the m — > oo limit. Explicit calculation [13] gives for i?^ = G{m/t^) 

m ^ ib ( Rl'\l-R^) V 
g^'P \{c{l - R,) - 4bR,Y ) 

with b, c defined as before and satisfy 46/c = 2 — e = 1/(3. Clearly, for large m, the denominator on the right 
hand side vanishes giving 

In the linear a-model eq.(7.6) is slightly modified. The masses are given by 

M'^ = -t+^v\ M', = -t+^v\ M'p=g\^ (7.9) 

D 2 

where t is the curvature of the scalar potential {t > corresponding to the broken phase) and v =< a >. 
The equation of state is h = vM^. One can obtain the relation between the mass ratios from eq.(7.9). The 
analogue of eq.(7.6) in this case is 
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-Rtt — 1 — 



A 



(7.10) 



12^2 



Here, the universal curve is a straight line again, but with different slope determined by the ratio of the 



four-fermi is a special case of the a model (as is well known) - for g"^ = 12A they describe the same low energy 
physics. Since both mass ratios are low energy quantities, the slope should be in fact a ratio of renormalized 
couplings (of course, this is not visible in the MF treatment). In other words, the curve i?,^ = Rt^{Rf) 
should have no knowledge of the bare parameters. In the four-fermi theories the renormalized couplings are 
independent of the bare four-fermi coupling G and the ratio is simply one. The reason behind M„ = 2Mp 
is the relative magnitude of A and g. From eq.(7.9) we see that the magnitude of A determines the size 
of the mass scale v (and meson masses), whereas g gives the magnitude of the fermion mass in units of v. 
In four-fermi theory, unlike the cr-model, mesons are composite rather then elementary. The compositeness 
condition forces the ratio of the couplings to be such that the cr is a true bound state i.e. < 2Mp. 
This requirement places a bound on the ratio g^ /X and is intimately related to the compositeness of the a 
and TT. If one looks closely at the two couplings A and g, two things become apparent. They appear in the 
lagrangian as gcyipt/j and Act'*. If the a is composite, then it renormalizes the same way as ipip. Therefore, 
compositeness requires that g^ and A renormalize the same way. So, if A^ = ZX, then = Zg^ and 
Xn/gji = A/(?2. Thus, radiative corrections cancel and the slope of the plot is independent of them. This, in 
fact, is what must happen if the scalars are composites because the wavefunction renormalization vanishes. If 
one accepts the democratic principle that, in an interacting theory, everybody has to interact with everybody 
else (unless there are selection rules that forbid it), then in order for both couplings to be non-zero, they 
have to renormalize the same way. In the cr-model scalars are pointlike and the ratio g^/X can assume any 
value. 

In the cr-model, the ratios are not constrained a priori since all the degrees of freedom are pointlike. 
Nevertheless, the largc-m behavior is consistent with that of the four-fermi model indicating that composit- 
ness is not the crucial ingredient here. Rather, it is chiral symmetry and the scaling of the mass ratios. In 
terms of the bare parameters, t, h, the mass ratio reads 



which, apart from the different critical exponents, is the same as eq.(7.7) obtained for the four-fermi theory. 
Prom here it follows that, in the /i — > oo limit, R„ ^ 1/3 = 1/S as it should. 

It has been argued many times in the literature that either the a model or the four-fermi model can be 
used to study the low energy regime of QCD. This is correct provided the bare masses are sufficiently small. 
These models are based on chiral symmetry and should be representative of QCD as long as the bound states 
that are studied are collective. Neither model is capable of describing atomic quarkonia like charmonium 
etc. In the heavy quark limit the predictions of these models should be qualitatively different from QCD. 
Whether or not they are suitable for the strange quark sector needs further investigation. It is interesting 
to note that an upper bound on the pseudo-Goldstone mass follows from the scaling plot eq.(7.10). In the 
broken phase, it is easy to see, that the following inequality holds 



coupling constants. For a fixed value of X/g"^ a different value of Rp emerges in the chiral limit. In that sense 




(7.11) 
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< -^^^F (7.12) 

In four dimension (5 = 3) gives M„ < \f2Mp. This bound is universal and in four dimensions holds to all 
orders in 1/A'' since the exponent 5 does not receive any corrections. The fact that an upper bound, m — > oo, 
on the TT — cr mass ratio is given by 1/^ < 1, shows that an explanation of this splitting as an effect of the 
standard spin-orbit interaction (a la atomic models) is problematic simply because the splitting survives 
even in the infinite mass limit. 



8. CONCLUSIONS 

In this paper we have illustrated how correlation length scaling and chiral symmetry together constrain 
and simplify the universal features of chiral symmetry breaking phase transitions. Perhaps the most useful 
results were 1. the explicit formula for in terms of the universal function /, 2. the relation between 
the pion mass, the chiral condensate, and the critical index 77, and 3. the simple properties of the universal 
function f itself. Within the context of the analysis of simulation data, each observation 1.-3. should lead to 
independent determinations of critical couplings and indices. Consistent results from all the methods should 
comprise convincing evidence that one has found universal features in the model of interest. 

We have chosen to illustrate these ideas in several simple models such as mean field theory and an 
ultra-violet fixed point with power-law singularities. Other examples can also be considered and the general 
approach of this paper can be applied. In ref.[22], for example, we considered both a fixed point and a 
logarithmically trivial sigma model to describe the data of four flavor lattice QED. Although both models were 
able to fit the equation of state, only the fixed point model gave a consistent description of the spectroscopy 
data for R.^. This exercise illustrates nicely how the use of PC AC and scaling in unison can help lead us to 
the correct theoretical scenario even when a comprehensive microscopic theory of a critical point is missing. 

We look forward to developing the ideas of this paper further and applying them to other presently 
mysterious chiral phase transitions such as QCD at nonzero chemical potential and four fermi models with 
continuous chiral symmetries and few fiavors in less than four dimensions. 
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FIGURE CAPTIONS 



1 . i?7r {t, m) as a function of m at fixed t in the symmetric and broken phase. The horizontal hne corresponds 
to the critical point f = 0. 

2. Typical behaviour of the universal ratio f{x)/f{0) as a function of {x/xq), x being the scaled variable 

t/ < v^v >^/^. 

3. The linear dependence of on < ■iptp >^ for a mean field theory. 

4. Expected behaviour of versus < ■^V >^ in a theory with anomalous dimension. 
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